Energy dependence of p/p ratio in p-\-p collisions 
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We have compiled the experimentally measured p/p ratio at midrapidity in p+p collisions from 
y/s = 23 to 7000 GeV and compared it to various mechanisms of baryon production as implemented 
in PYTHIA, PHOJET and HIJING/B-B models. For the models studied with default settings, 
PHOJET has the best agreement with the measurements, PYTHIA gives a higher value for y/s < 
200 GeV and the ratios from HIJING/B-B are consistently lower for all the y/s studied. Comparison 
of the data to different mechanisms of baryon production as implemented in PYTHIA shows that 
through a suitable tuning of the suppression of diquark-antidiquark pair production in the color 
field relative to quark-antiquark production and allowing the diquarks to split according to the 
popcorn scheme gives a fairly reasonable description of the measured p/p ratio for y/s < 200 GeV. 
Comparison of the beam energy dependence of the p/p ratio in p+p and nucleus-nucleus (A+A) 
collisions at midrapidity shows that the baryon production is significantly more for A+A collisions 
relative to p+p collisions for y/s < 200 GeV. We also carry out a phenomenological fit to the t/beam 
dependence of the p/p ratio. 
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Protons (j>) and anti-protons (p) are the most 
abundantly produced baryons in high energy colli- 
sions. These have been measured at various cen- 
ter of mass energies (y/s) in hadron-hadron [J 0] 
and nucleus-nucleus collisions [3j as a function of ra- 
pidity (y) and transverse momentum (pr). Rapid- 
ity dependence of baryon production is expected to 
provide information on baryon transport and stop- 
ping Q and the pt dependence of the yields is ex- 
pected to help in understanding the baryon pro- 
duction mechanism |2|. The p/p ratio within the 
assumption of a thermal model is used to obtain 
the baryon chemical potential in heavy-ion colli- 
sions 5, 6]. Recently it has been argued based on 
QCD that there are constraints on allowing quarks 
to trace the baryon number Q , although they carry 
a baryon number of 1/3 based on quark model clas- 
sification. It is argued that the trace of the baryon 
number could be associated with non-perturbative 
configurations of gluon fields rather than to valence 
quarks. All these make the study of the energy de- 
pendence of p/p ratio important in high energy col- 
lisions, where the role of the gluonic contributions 
to particle production is expected to increase with 

In the string picture the process of baryon pro- 
duction is not unique. Mesons in such a picture can 
be viewed with a short string between a quark and 
anti-quark endpoints. However for the baryons con- 
sisting of three quarks it is difficult to visualize in a 
simple way. Baryon production in string picture is 
implemented in the PYTHIA model |8j. The sim- 
plest mechanism of baryon production in such a pic- 
ture is through a diquark model. Any quark of a 
given flavor is assumed to be represented either by a 



quark or an antidiquark in a color triplet state. Then 
the baryon and antibaryon are produced as nearest 
neighbors along the string. Such a model has to deal 
with the relative probability to pick a diquark over a 
quark. The extra suppression associated with a di- 
quark containing a strange quark purely from phase 
space considerations and when a baryon is formed by 
joining a diquark and a quark, it has to be a sym- 
metric three-quark state. Another equivalent mech- 
anism is that in which diquarks as such are never 
produced, but baryons appear due to the successive 
production of several quark-antiquark pairs. Such a 
mechanism is referred to as the popcorn mechanism. 
These pairs exist by means of the color fluctuations 
in the field [9|. An advanced version of the pop- 
corn mechanism is described in Ref . [l(} ■ While the 
simpler popcorn mechanism admits at most one in- 
termediate meson formation, the advanced version, 
on the other hand, allows for the possibility of many 
such mesons. 

Another model of particle production which is 
widely used for comparison to data is the PHO- 
JET 53. The PHOJET, is a two component model 
that combines the ideas of the Dual Parton Model 
(DPM) [l2| (soft processes) with perturbative QCD 
(hard processes). The mechanism of Pomeron ex- 
change is at the heart of the DPM. According to 
the DPM, the leading contribution to multiparticle 
production in high-energy hadron-hadron collisions 
arises from the exchange of a single Pomeron be- 
tween the colliding hadrons. Secondary Pomeron 
exchanges account for the remaining activity in the 
event. Each exchanged Pomeron gives rise to two 
color-neutral chains stretching between quarks and 
diquarks, for baryons, or quarks and anti-quarks for 
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mesons. The basic difference between PYTHIA and 
PHOJET lies in their approach towards an event 
formation. The starting point of particle production 
in PYTHIA is through the description of possible 
hard interactions in e + + e _ , p+p{p) or e+p col- 
liders and then combines several ideas for the soft 
hadronic interactions, whereas in PHOJET model 
it initializes the event generation by describing the 
soft component of hadron-hadron, photon-hadron or 
photon-photon interactions at high energies. The 
hard component is introduced later and calculated 
by perturbative QCD at the partonic level. 

A novel mechanism of baryon transport motivated 
by the Regge theory [l3| and differing from the di- 
quark breaking model has been implemented in the 
form of an event generator, the HIJING/B-B [Til ]. 
As discussed briefly earlier, the mechanism is moti- 
vated from the non-perturbative gluon field config- 
uration called the the baryon junction. The baryon 
junction is found to be originating from the basic 
concepts of QCD and is a vertex where the color 
flux lines flowing from the three valence quarks are 
connected. The junction is expected to play a dy- 
namical role through the Regge exchange of junc- 
tion states in high energy collisions. The junction 
exchange could provide a natural mechanism for the 
transport of baryon number into the central rapidity 
region. Further details can be found in Ref. @, [lij 

In this article we compare the energy dependence 
of the experimentally measured p/p ratio at midra- 
pidity in p+p collisions for various to the above 
discussed models. The data for the i/s = 23, 31, 45, 
53 and 63 GeV are from the ISR experiments 
The p/p ratio at ^/s = 200 GeV is from the STAR 
experiment at RHIC @. And the p/p ratio at the 
highest energies of 900 GeV and 7 TeV are from AL- 
ICE at LHC @. 

Figure [1] shows the p/p ratio at midrapidity for 
various yfs (except y/s = 63 GeV, reasons for not 
showing these are discussed later), in p+p collisions 
as a function of Pt- All data points are only shown 
for the pt < 1 GeV/c. We observe that the ratios 
are constant as a function of pt for each t/s and the 
value of the p/p ratio increases with ^/s. The solid 
lines show the recent measurements by the ALICE 
experiment at LHC [2j. The Iowbt dependence of 
the p/p ratio is discussed in Ref. [2}, where a drop of 
p/p with increase in pr is expected in HIJING/B-B. 
In the this paper we discuss only the energy depen- 
dence of the p/p ratio at midrapidity by comparison 
to various models of baryon production. 

We first compare the experimental measurements 
to PYTHIA, PHOJET and HIJING/B-B models, 
then compare the p/p ratio to various baryon pro- 
duction mechanism as implemented in PYTHIA. 
This is followed by the discussion of asymmetry in 
proton and anti-proton production and comparison 
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FIG. 1: p/p ratio at midrapidity as a function of pt for 
p+p collisions. The solid lines corresponds to new data 
from ALICE at LHC @] . The p/p ratio for y/s = 200 GeV 
is from the STAR experiment at RHIC |{J. The data for 
V^s < 200 GeV are from the ISR [J]. The dashed lines are 
straight line fits to the px dependence of the p/p ratios 
at various y/s, assuming the ratios do not depend on pt 
for the measured range. 



to corresponding available results from heavy-ion 
collisions. 

Figure [2] shows the increase of the pt integrated 
p/p ratio at midrapidity with increase in ^fs for p+p 
collisions. The experimental data are compared to 
three models,viz, PYTHIA (Ver. 6.4), PHOJET 
(Ver. 1.12) and HIJING/B-B (Ver. 1.34),all with 
default settings. The p/p ratio for y/s — 63 GeV 
shows an abnormally high value, although shown in 
the figure, we would not consider it for physics dis- 
cussions. It is expected that RHIC data collected in 
the year 2005 at y/s = 63 GeV will help in resolving 
the abnormality in the p/p ratio. All models studied 
show that the p/p ratio increases with and ap- 
proaches unity for higher energies (LHC). Infact, the 
PYTHIA and PHOJET models give very similar val- 
ues at the LHC energies, while HIJING/B-B under 
predicts the p/p ratio. The major difference occurs 
for y/s < 200 GeV, PYTHIA model gives higher val- 
ues and HIJING/B-B continues to give lower values 
of p/p ratio. Only the PHOJET model with de- 
fault settings gives a reasonable description of the 
y/s dependence of the measured p/p ratio for p+p 
collisions. 

The PYTHIA model has some variations in the 
baryon production mechanism. In Fig. [3] we com- 
pare the experimental p/p ratio to such variations as 
implemented in the model. The parameter that we 
varied is known as MSTJ(12), it can take up values 
from to 5, with the value of 2 as the default setting 
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FIG. 2: (Color online) p/p ratio at midrapidity as a func- 
tion of y/s for p+p collisions. The experimental data are 
compared to model calculations from PYTHIA [§] , PHO- 
JET [n[ and HIJING/B-B [jj] with default settings. 
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FIG. 3: (Color online) p/p ratio at midrapidity as a 
function of y/s for p+p collisions compared to various 
implementation of the baryon production schemes in 
PYTHIA. See text for more details. 
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FIG. 4: (Color online) Asymmetry for proton and anti- 
proton production at midrapidity for p+p collisions as a 
function of y/s. The solid line is a fit to the data, with 
the functional form shown. 



(one used in FigHJ}. We did not consider MSTJ(12) 
= 0, as it corresponds to no baryon-antibaryon pair 
production. The condition MSTJ(12) = 1 refers to 
the mechanism where baryon production is through 
diquark-antidiquark pair production with the di- 
quark being treated as a unit. While MSTJ(12) = 2 
has the additional possibility for diquark to be split 
according to the popcorn scheme. The mechanism 
of baryon production for the case MSTJ(12) = 3 
is same as that for MSTJ(12) = 2, but has an ad- 
ditional condition that the production of first rank 
baryons may be suppressed. For this case, we addi- 
tionally changed the value of the parameter which 
governs the suppression of diquark-antidiquark pair 
production in the color field, compared with quark- 
antiquark production. The value we put in for this 
parameter is 0.05 compared to the default value of 
0.1. This is referred to as MSTJ(12)=3-Tuned in 
the Fig. 13 The condition MSTJ(12) = 4 again re- 
volves around MSTJ(12) = 2 with an extra condition 
that the diquark vertices are suppressed. The last 
scheme implemented corresponds to MSTJ(12)= 5 
is similar to MSTJ(12) = 2, but with an advanced 
version of the popcorn model. Fig. [3] shows that 
for the lower beam energies the model results with 
the condition MSTJ(12)=3- Tuned has a reasonable 
agreement with the experimental data. For higher 
y/s all the above conditions give similar values oip/p 
ratio. 

Different baryon production mechanisms could 
lead to an asymmetry in the production of protons 
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FIG. 5: (Color online) Comparison of beam energy of p/p 
ratio at midrapidity for p+p and nucleus-nucleus (A+A) 
collisions. The solid, dash and dotted lines are fits to 
the p+p and A+A data inspired by the model based on 
baryon string junction picture 0, El- 



and anti-protons. This asymmetry can be measured 
by constructing the following ratio, 



N p 



Np + Np' 



(1) 



where N p and Np are the number of protons and 
anti-protons. As pair production would lead to same 
number of protons and anti-protons, the asymmetry 
will have a value of zero. Any non-zero value indi- 
cates the fraction of protons in midrapidity due to ef- 
fects such as stopping. Figure |4] shows the asymme- 
try ratio for protons and anti-protons as measured in 
p+p collisions for various y/s ranging from 23 GeV to 
7 TeV. The asymmetry is found to decrease with in- 
crease in y/s, indicating the decreasing contributions 
of protons due to stopping at midrapidity. The ratio 
changes from about 46% at y/s = 23 GeV to 0.5% at 
the top LHC energy of 7 TeV. This range in y/s cor- 
responds to a range in ?/boam of 3 to 9 units, respec- 
tively. This information is useful to study double 
baryon production in p+p collisions Q and baryon 
number flow over long rapidity interval (H)]. The 
solid line in Fig. [4] is a fit (x 2 /ndf = 3/4) to the ex- 
perimental data with the function Ae~ BL ° 9 ^ , with 
the parameters A = 6.7 ± 0.9 and B = 0.85 ± 0.04. 

In Fig. [5] we compare the beam energy depen- 
dence of the experimentally measured midrapidity 
p/p ratio in p+p collisions to the available p/p ratios 
at midrapidity in nucleus-nucleus (A+A) collisions. 
The A+A collision data are taken from the exper- 
iments at AGS 0, SPS [13 and RHIC [3]. For 
both the systems the p/p rapidly rises with beam 



energy and approaches unity. For the p+p collisions 
the p/p ratio has a value of 0.991 ± 0.005(stat.) ± 
0.014(syst.) at 7 TeV, while for heavy-ion collisions 
the p/p ratio has a value of 0.77 ± 0.14 at 200 GeV. 
Looking at the region in beam energy where there 
is overlap between p+p and A+A collisions, the rel- 
ative proton contributions at midrapidity for A+A 
collisions is more than for p+p collisions. The val- 
ues of the ratio seem to become equal around 200 
GeV. Also shown in the figure are fits to the exper- 
imental p/p ratio in the p+p collisions to a func- 
tion of the form [1 + Cexp[(aj — ap)?/bcam] 1 • The 
dashed line corresponds to aj — ap = (1.2 - 0.5) = 
-0.7, values as expected from a Regge model where 
the baryon pair production at very high energy is 
governed by Pomeron exchange and baryon trans- 
port by string-junction exchange 0, fl8| . The aj 
and ap parameters corresponds to the junction in- 
tercept and the Pomeron intercept in the models. 
The x 2 /ndf for the fit is 66/6. The best fit (x 2 /ndf 
= 4/5) is however obtained for aj — ap — -0.97 ± 
0.05 (solid line), with about a factor of 3 increase 
in the value of parameter C compared to the results 
presented in dashed lines. The value of C for the 
case is 10.3 ± 0.4 and that for the second case (solid 
line) is 35 ± 7. The Regge model inspired fit to the 
heavy-ion data (dotted line), yields a poor \ 2 l n df = 
33/7, with fit parameters C = 119208 ± 49600 and 
a j — ap = -2.9 ± 0.2. The fit misses the the higher 
energy data points where such a model is more reli- 
able. Constraining the fit to the energy range of 23 
to 200 GeV gives C = 205 ± 580 and aj - a P = 
-1.2 ± 0.6 with \ 2 /ndf — 0.6/1. Upcoming heavy- 
ion data at higher beam energies in LHC will give 
a clear picture of applicability of such a model to 
heavy-ion collisions. 

In summary, we have presented a compilation of 
the available data for p/p ratio at midrapidity for 
p+p collisions as a function of y/s. We have also 
compared these ratios to the beam energy depen- 
dence from heavy-ion collisions and found that be- 
low y/s~tm = 200 GeV, the proton contribution at 
midrapidity in A+A collisions is significantly more 
compared to those in p+p collisions. The p/p ratio is 
constant as a function of px for all the beam energies 
forpT < 1 GeV/c. This experimental observation al- 
ready puts a constrain on mechanism of baryon pro- 
duction such as those implemented in HIJING/B-B. 
We also compared the p/p ratio vs. y/s to results 
from various models with different baryon produc- 
tion mechanisms, such as PYTHIA, PHOJET and 
HIJING/B-B with default settings. It is observed 
that PHOJET gives the best description of the data 
for all y/s, PYTHIA gives higher values of the p/p 
ratio for y/s < 200 GeV and HIJING/B-B under 
predicts the ratio for all beam energies. A detailed 
investigation of various mechanisms of baryon pro- 
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duction as implemented in PYTHIA shows that the 
baryon production through diquark-antidiquark pair 
production with the diquark being treated as a unit 
and the additional possibilities (arrived by tuning 
various parameters) of diquark splitting according 
to the popcorn scheme and the production of first 
rank baryons suppressed gives a reasonable descrip- 
tion for the p/p ratio for yfs < 200 GeV. The asym- 
metry a measure of proton stopping at the midra- 
pidity in p+p collisions are presented. The frac- 
tion of protons stopped around midrapidity varies 
from 46% at yfs = 23 GeV to 0.05% for = 7 
TeV. This energy range corresponds to a range in 
2/bcam from 3 to 9 units, respectively. The data has 
also been compared to baryon string junction moti- 



vated phenomenological function whose parameters 
can constrain the Regge-model inspired descriptions 
of baryon asymmetry in p+p collisions. 
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